Introduction

Simulate some mother-offspring data and illustrate partitioning of variance into additive and environmental components using rstanarm assuming that offspring are half siblings

plot(cars)

Lab 12

As part of her dissertation, Nora Mitchell measured

Protea punctata

Protea venusta

  1. What is the heritability of each trait in each species? (3 traits, 2 species)

  2. What is the estimated selection differential in each species? (3 traits, 2 species)

  3. What is the predicted response to selection in each species? (3 traits, 2 species)

Just for me while I’m getting this ready

library(tidyverse)

rm(list = ls())

field <- read_csv("Protea_field.csv")
Rows: 147 Columns: 5
── Column specification ─────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): species
dbl (4): seedheads, lma, fwc, lwr

ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
table(field$species)

punctata  venusta 
      86       61 
greenhouse <- read_csv("Protea_greenhouse.csv")
Rows: 437 Columns: 5
── Column specification ─────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): species
dbl (4): mom_id, lma, fwc, lwr

ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
table(greenhouse$species)

punctata  venusta 
     245      192 
LS0tCnRpdGxlOiAiRXN0aW1hdGluZyBoZXJpdGFiaWxpdHkgZnJvbSBoYWxmLXNpYmxpbmdzIgpvdXRwdXQ6IAogIGh0bWxfbm90ZWJvb2s6CiAgICB0b2M6IHllcwogICAgdG9jX2Zsb2F0OiB0cnVlCi0tLQoKIyBJbnRyb2R1Y3Rpb24KClNpbXVsYXRlIHNvbWUgbW90aGVyLW9mZnNwcmluZyBkYXRhIGFuZCBpbGx1c3RyYXRlIHBhcnRpdGlvbmluZyBvZiB2YXJpYW5jZSBpbnRvIGFkZGl0aXZlIGFuZCBlbnZpcm9ubWVudGFsIGNvbXBvbmVudHMgdXNpbmcgYHJzdGFuYXJtYCBhc3N1bWluZyB0aGF0IG9mZnNwcmluZyBhcmUgaGFsZiBzaWJsaW5ncwoKYGBge3J9CnBsb3QoY2FycykKYGBgCgojIExhYiAxMgoKQXMgcGFydCBvZiBoZXIgW2Rpc3NlcnRhdGlvbl0oaHR0cHM6Ly9vcGVuY29tbW9ucy51Y29ubi5lZHUvZGlzc2VydGF0aW9ucy8xMzcwLyksIE5vcmEgTWl0Y2hlbGwgbWVhc3VyZWQKCiFbX1Byb3RlYSBwdW5jdGF0YV9dKFByb3RlYV9wdW5jdGF0YS5wbmcpCgohW19Qcm90ZWEgdmVudXN0YV9dKFByb3RlYV92ZW51c3RhLnBuZykKCgoxLiBXaGF0IGlzIHRoZSBoZXJpdGFiaWxpdHkgb2YgZWFjaCB0cmFpdCBpbiBlYWNoIHNwZWNpZXM/ICgzIHRyYWl0cywgMiBzcGVjaWVzKQoKMi4gV2hhdCBpcyB0aGUgZXN0aW1hdGVkIHNlbGVjdGlvbiBkaWZmZXJlbnRpYWwgaW4gZWFjaCBzcGVjaWVzPyAoMyB0cmFpdHMsIDIgc3BlY2llcykKCjMuIFdoYXQgaXMgdGhlIHByZWRpY3RlZCByZXNwb25zZSB0byBzZWxlY3Rpb24gaW4gZWFjaCBzcGVjaWVzPyAoMyB0cmFpdHMsIDIgc3BlY2llcykKCgojIyBKdXN0IGZvciBtZSB3aGlsZSBJJ20gZ2V0dGluZyB0aGlzIHJlYWR5CgpgYGB7cn0KbGlicmFyeSh0aWR5dmVyc2UpCgpybShsaXN0ID0gbHMoKSkKCmZpZWxkIDwtIHJlYWRfY3N2KCJQcm90ZWFfZmllbGQuY3N2IikKdGFibGUoZmllbGQkc3BlY2llcykKZ3JlZW5ob3VzZSA8LSByZWFkX2NzdigiUHJvdGVhX2dyZWVuaG91c2UuY3N2IikKdGFibGUoZ3JlZW5ob3VzZSRzcGVjaWVzKQpgYGA=